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1. Introduction 



Let fn{t) be a real random polynomial of degree n, 

/n(t) = Co + Cit+--- + Cnr, (1.1) 

where cq, ci, . . . , are independent real random variables. Distribution of zeros for various 
classes of random polynomials is studied in the classical papers by Bloch and Polya [BP], 
Littlewood and Offord [LO], Erdos and Offord [EO], Erdos and Turan [ET], and Kac 
[K1-K3]. We will assume that the coefficients cq, ci, . . . , c„ are normally distributed with 

Ecj = 0, Ec2 = (7|. (1.2) 

In the case when 

fn{t) is the Kac random polynomial. Another interesting case is when 

As is pointed out by Edelman and Kostlan [EK], "this particular random polynomial is 
probably the more natural definition of a random polynomial". We call this polynomial 
the SO (2) random polynomial because its m-point joint probability distribution of zeros 
is S0(2)-invariant for all m (see section 5 below). The SO (2) random polynomial can 
be viewed as the Majorana spin state [Maj] with real random coefficients, and it models 
a chaotic spin wavefunction in the Majorana representation. See the papers by Leboeuf 
[Lebl, Leb2], Leboeuf and Shukla [LS], Bogomolny, Bohigas, and Leboeuf [BBL2], and 
Hannay [Han], where the SU(2) and some other random polynomials are introduced and 
studied, that represent the Majorana spin states with complex random coefficients. 

Let {ti, . . . , Tfe} be the set of real zeros of fn{t)- Consider the distribution function of 
the real zeros, 

Pnit) = Ei^{j:Tj<t}, 

where the mathematical expectation is taken with respect to the joint distribution of the 
coefficients cq, . . . , c^. Let 

p^(t) = p;(t) 

be the density function. By the Kac formula (see, e.g., [K3]), 



_ V^.(t)C.(t) - Bljt) 
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where 



n 



Kit) 



(1.4) 



i=i 

n 



The derivation of (1.3) by Kac is rather complex. A short proof of (1.3) is given in the 
paper [EK] by Edelman and Kostlan. See also the papers by Hannay [Han] and Mesincescu, 
Bessis, Fournier, Mantica, and Aaron [M-A], and section 2 below. The formula (1.3) implies 
that for the Kac random polynomial, 



lim pn{t) =p{t) = 



1 



and 



Pn(±l) = 



7r|l-t2|' 
n(n + 2) 



(1.5) 



12 



1/2 



(see [K3], [BS], and [EK]). The limiting density p(t) is not integrable at ±1, and this means 
that the zeros are mostly located near ±1. Observe, in addition, that Pn{t) is an even 
function of t, and the distribution pji{t)dt is invariant with respect to the transformation 
t — > 1/t. Kac [Kl] proves that the expected number of real zeros has the asymptotics 

f oo 



/oo 
Pn{t)dt={2/7r)\ogn + 0{l). 
-oo 



Kac [K2], Erdos and Offord [EO], Stevens [Ste], Ibragimov and Maslova [IM], Logan and 
Shepp [LS], Edelman and Kostlan [EK], and others extend this asymptotics to various 
classes of the random coefficients {cj}. Maslova [Masl] evaluates the variance of the 
number of real zeros as 



Var#{j:/„(r,)=0} 



TT 



1-- 1 lnn(l + o(l)), 

TT 



n 



oo, 



and she proves the central limit theorem for the number of real zeros (see [Mas2]), for a 
class of distributions of the random coefficients {cj}. 

In this paper we are interested in correlations between the zeros Tj of the Kac random 
polynomial. Let us consider first the zeros in the interval (—1, 1). Define straightening of 
Tj as 



= P(t,), P{t) 



p{u) du. 



3 



In the limit when n — > oo, the straightened zeros Q are uniformly distributed on the real 
line, so that 

lim Ei^{j: a<Cj <b} = b-a. (1.6) 



Prom (1.5) we get that 



P(t) = / — — = — In 

^ ^ Jo TT 1 - «2 27r 



1 + t 



1-t 



1 

= — artanh t, 

TT 



hence ^ 

= — artanh Tj . (l-^) 

Let Pmn{si, . . . , Sm) be the joint probability distribution density of the straightened zeros 

PnmySi, . . . , Sm) — iini I a a I 

Asi,...,As^^O \Asi...ASm\ 

(1.8) 

It coincides with the correlation function 

, . s E [^ri{si,Si+Asi)...CniSni:Sm + ^Sm)] 

knm{Si,...,Sm)= ^ lim — ■ . (1.9) 

Asi,...,As^^O \Asi...ASm\ 

where 

Cn{a,b) = #{j:a<Q<b}. 

We assume in (1.8) and (1.9) that Si ^ Sj for all i ^ j. Our aim is to find the limit 
correlation functions 

km{si,...,Sm) ^ lira knm{si,...,Sm)- (1-10) 
n—^oo 

We prove the following results. 

Theorem 1.1. The limit two-point correlation function k2{si,S2) of the straightened 
zeros (j — artanh Tj- of the Kac random polynomial is equal to 

7 / \ i2 / \ I sinh7r(si — 82)1 . 1 

«;2(si,S2) = tanh nlsi - S2) ^ 5 — ; arcsm — - — (l-H) 

cosh 7r(si - S2) cosh7r(si - S2) 



Observe that /c2(si, S2) depends only on si — S2, and it has the following asymptotics: 

^2(^1, S2) = —\si- S2I + 0{\si - S2p), |si - S2I 0, 

k2{su S2) = 1 - ^ e-^^^l^i-^^l + 0(e-^'^l^i-*^l), |si - S2I ^ 00. 



1- 



0.6- 



0.4- 



0.2- 




Fig 1: The two-point correlation function of straightened zeros of the Kac random polynomial. 

The graph of ^2(0, s) is given in Fig. 1. 

Theorem 1.2. The limit m-point correlation function km{si, . . . , Sm) of the straight- 
ened zeros — TT~^ artanhTj is equal to 



/oo l-OO 
■■I \yi...ym\e~^^^'^"''^^dyi...dym, 
-00 J —00 



where Y = (j/i, . . . , ym) and the matrix Vm is defined as 

1 



(1.12) 



(1.13) 



In particular, km{si, . . . , Sm) depends only on the differences of si, . . . , Sm, hence it is 
translation invariant. 

The proof of Theorems 1.1 and 1.2 is given in sections 2, 3 and 4 below. It is based 
on computation of the determinant of some matrices which consist of 2 x 2 blocks. This 
computation is of independent interest. The basic example is the matrix 



where 



1 



/All A12 
A21 A22 

V A^i A^2 



1 titj (1 ^i^j) 
tj 1 I tj^tj 



Ai^ \ 

A2m 



(1.14) 



, m. 



(1.15) 



\ (1 - titj)^ (1 - titj)^ 
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We prove in section 4 that 

det = 



<i<j<m 



{ti tj 



(1.16) 



It is interesting to note that determinants of matrices consisting of 2 x 2 blocks appear 
also in the theory of random matrices (see, e.g., [Dys] and [Meh]), statistical physics, and 
other applications. 

Consider now zeros tj with |tj| > 1. Define straightening of tj as 



p{u) du if t < —1, 

= P(r,), P{t)^{ 

it ^'^^^ 



(1.17) 



du if t> 1. 



In the limit when n — cxd, the straightened zeros Q are uniformly distributed on the real 
line. From (1.5) 



1+t 



1-t 



1 



— artanh t ^ , 



TT 



SO that 



G = TT ^ artanh T- ^ . 



(1.18) 



(1.19) 



Denote by A;°^(si, . . . , Sm) the correlation function of the straightened zeros Qj with |tj | > 
1. 



Theorem 1.3. 



^nm (^1 ) • • • ) ^m) — ^nm(si ) • • • ) ^m) • 



(1.20) 



In other words, the correlation functions of the straightened zeros outside of the 
interval (—1, 1) coincide with those inside of the interval. Finally, let us consider correlation 
between zeros inside of the interval (—1, 1) and outside of this interval. Let Knm{tii ■ ■ ■ i tm) 
be the correlation function of the zeros Tj (without straightening). 

Theorem 1.4. Assume that \ti\, . . . ,\ti\ < 1 and . . . , \tm\ > 1- Then the limit 

lim Knmih, ■ ■ ■ ,tm) = Km{ti, . . . ,tm) (1-21) 
n— >oo 

exists and 

Km{tl, • • • , tm) = Ki{ti, .. .,ti) Km-l{tl+i, . . .tm)- (1.22) 
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This means that the zeros inside and outside of the interval (—1,1) are asymptotically- 
independent. Observe that 

(^1 5 • • • 5 

provided that either all \tj\ < 1 or all \tj\ > 1 (cf. the formula (2.14) below). Proof of 
Theorems 1.3 and 1.4 is given in the end of section 4. 

In sections 5 and 6 we investigate correlation functions of real zeros of the SO (2) 
random polynomial. 

2. General Formulae 



■^m (tl : • • • : ^m) 
p{ti) ...p{tm) _ 



(1.23) 



Let 

n 

/n(0 = E^^-^'' (2.1) 

be a polynomial whose coeflBcients Cj are random variables with an absolutely continuous 
joint distribution. Let 

^n{a, h) = #{Tfc :a<Tk<b, U{Tk) = 0} (2.2) 

be the number of real roots of fn{t) between a and 6, and let Pn{t) be the density of real 
zeros tk of fn{t), so that 

ECn{a,b)^ [ pn{t)dt. (2.3) 

J a 

It is not difficult to show that 

/oo 
\y\DniO,y;t)dy, (2.4) 
-oo 

where Dn{x,y;t) is a joint distribution density of /n(t) and fn{t), 

nh nd 

VT{a<fr,{t)<h;c<f'^{t)<d}^ / Dn{x,y-t)dxdy. (2.5) 

J a J c 

Indeed, if f'n{t) = y then asymptotically as At — > 0, the function fn{t) has a zero in the 
interval [t,t + At] if fn{t) is in the interval [0,—yAt], and this gives (2.5). Similarly, the 
m-point correlation function Knmiti, ■ ■ ■ ,tm) for pairwise different ti, . . . ,tm is equal to 

/oo poo 
•■■ / \yi ■ ■ - ymlDnmi^jyi, ■ ■ ■ ,Q,ym'iti, . . . ,tm)dyi . . .dym, 
-oo J —oo 

(2.6) 



where Dnmi^ijUij ■ • • j ^miym'i til ■ ■ ■ 1 tm) is a joint distribution density of the vector 

Fn = {fn{tl), fni^l), ■ ■■, fn{tm),fn{tm)), 

SO that 

Pr {ai < fn{tl) < h; d < fn{ti) < di; . . . ; am < fn{tm) < hm\ Cm < fn{im) < dm} 

"bl pdl />bm fdm 

■Dnm i-^l: VIt ■ ■ ■ T '^rm Urm tlj ■ ■ ■ j tm) dX\dy\ . . . dXm,dym- 



(2.7) 

If {cj} are independent random variables with Var cj > then the covariance matrix 
of the vector is positive, provided that n > 2m — 1 (see Appendix B at the end of 
the paper). Similar formulae are derived for the correlation functions of complex zeros of 
random polynomials with complex and real coeflBcients (see [Han] and [M-A]). 

Observe that 

™- rbi pbm 

^ Yi^n{aj,bj) = •••/ Knm{tl,---,tm)dti...dtm, (2.8) 

j — l J 0,\ ''dm 

provided that (ai, 61), . . . , (a^, hm) are pairwise disjoint, and 

Pnit) = Kr^iit), E(en(a,6))= / Kniit)dt. (2.9) 

J a 

For the general case, when (ai, bi), . . . , {am, bm) may intersect, we have the following ex- 
tension of (2.8): 



I 



( \ 



^Win{a,,b,)^ J2 n 

J = l {Ai,...,Ai)j = l 



dtj 



n {ai,bi) 



Knl{ti,...,ti), (2.10) 



where the sum is taken over all possible partitions (Ai, . . . , A^) of {1, . . . , m}, such that 

Ai nAj = ^, j, 

Ai[J---UAi = {l,...,m}, (2.11) 
\Ai\ > 1 1 = 1,. ..,1. 

In particular, when m = 2 we have 

E [Cn(oi,&l )Cn (02,62)] = / / Kn2{tlM)dtldt2 (2.12) 



if (ai, 61) n (02, 62) = 0, and 



E[enK&)]= / Pn{t)dt + 



-b 



Kn2{tl,t2)dt\dt2 



(2.13) 



From definition (1.9) of the m-point correlation function, it follows that the m-point cor- 
relation function knm{si, ■ ■ ■ iSm) of the straightened zeros Q = P{rj) is related to the 
m-point correlation function K^miti, ■ ■ ■ ,tm) of the zeros Tj by the formula 



^nm(^l) • • • ) ^m) — 



■^nm (^1 



5 ; 



p{tl) . . .p{tr 



(2.14) 



tl=P-l(si),...,t^=P-l(Sm) 



Assume now that the coefficients Cj are independent Gaussian variables with zero mean 
and the variances cr|, j = 0,...,n. Then Dni{x,y;t) is a Gaussian distribution density 
with the covariance matrix 



A = 



E fnit)m E (my 



2 - 



Bn{t) Cn{t) 



(2.15) 



where An{t), Bnit) and C^(t) are defined in (1.4), and from (2.4) we get the Kac formula 
(1.3). 



3. Two-Point Correlation Function for the Kac Polynomial 



Let fn{t) = Co + cit + ■ ■ ■ + Cnf^ be the Kac polynomial, so that Ck,k = 0, , 
real independent Gaussian random variables with 



, n, are 



Ecfc = 0, Ecfc2 = l. 



(3.1) 



Consider the covariance matrix A„ of the Gaussian vector (/n(^i), /n(*i)) /n(*2), /n(*2))- 
From (3.1) 

. _ 1 - iht2r+' 



E fn{tl)fn{t2) = J2(ht2y 



k=0 

E /;(tl)/n(^2) = ^ 

E /;(ti)/;(t2) = 



1 - tlt2 



n+1 



1 - {tlt2) 

1 - ht2 

1 - (tlt2)"+^ 



1 - ht2 



Assume that \ti\, \t2\ < 1. Then from (3.2) we obtain that 



(3.2) 



lim = A, 



(3.3) 



with - . 1 ^ 

Ji 1 ti \ 

. y^-ti^y l-tlt2 (l-tlt2)2 \ 

ti tz l+tita 

A = 



1 




(l-tl2)2 


1- 


-tlt2 


(l-tlt2)2 












l+tlt2 


(1- 


1 


l+tlt2 


(1- 


1 


(l-tlt2)^ 
t2 


1- 


-tlt2 


(l-tlt2)^ 


1 


-t2^ 


(l-t2^)^ 




tl 


l+tlt2 




t2 


l+t2^ 


(1- 


-tlt2)^ 


(l-tlt2)^ 


(1- 


-t2'y' 


(l-t2^)« 



(3.4) 



We prove in the section 4 below that 

det A = o,,/^'~^£' 7^ • (3.5) 

Let Q. be the two-by-two matrix obtained by removing the first and the third rows and 
columns from A~^. Then 

where 



A = 


(1 


-tit2)\l 








B = 


(1 








(3.7) 


C = 


(1 






-t2f 





By (2.6), the correlation function -^'2(^1,^2) is equal to 

-1 l-OO l-OO 

K2(tut2)= ^ , r-—- / / |yiy2|e-^(^"'^)dyidy2 (3.8) 
47r''vdet A J-00 J -00 

where Y = (yi, 2/2)- Since 

/ |2/i2/2|e-^(^^?+^^^^^^+^^^^)d|/idy2 
-00 ^ —00 



1 H , arcsm , = 



(3.9) 



^c(i - 52) V J ' 

(see Appendix A), we obtain that 

{ti-t2r 



K2{tlM] 



"7r2(l-M2)2(l-ti')(l-t2') 

^ Itl-t.l . v/(l-ti2)(l-t2^) 

H ■ — , arcsin -^^-^ — 

7r2(l - tit2)V(l-^l')(l-*2') 1 - ^1^2 



(3.10) 



Consider the correlation function A;2(-si, S2) of the straightened zeroes Qj — ^ artanhrj. 
By (2.14), 

A;2(si,S2) = /!^! 5 *i = tanh(7rsi), t2 = tanh(7rs2)- (3.12) 
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Since 

(see (1-5)), we obtain that 



= JT^^ + (TTT^^^? ""^^'^ 

= tanh^ 7r(si — S2) + ^ ' " ^ ^ arcsin 



cosh^7r(si - S2) cosh7r(si - S2) 

Theorem 1.1 is proved. 



4. Higher order correlation functions for the Kac polynomial 



Let fn{t) be the Kac polynomial, and let ti,t2, ■ ■ ■ ,tm he m > 3 distinct points in the 
interval (—1,1). Denote by A^-* the covariance matrix of the Gaussian vector 

(/n(tl),/;(tl),...,/n(tm),/;(tm)), 

and by the limit of A^^ as n — > 00, 

A^ = hm A^) (4.1) 



Then 



where 



and by (3.2), 



where 



[cf. (3.4)]. 



m 



,(n) 
ij 



11 
^21 



^12 
^22 



\A 



ml m2 



in) 



A H 



^fL{u)fn{tj) E/;(t,)/;(t,) 



/ An 
A21 



\A 



ml 



A12 
A22 

Am2 



Ai™ \ 
A2m 



1 ^iij (1 t-itjY^ 



(4.2) 



(4.3) 



(4.4) 



(4.5) 
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If iljn denotes the m x m matrix obtained by removing all the odd number rows and 
columns from A~^, then by (2.6), the correlation function Km,{ti, ■ ■ ■ , tm) is equal to 

Kmih, . . . ,tm) = -——j=== ■■■/ \yi...ym\e-^^^^^^^Uyi...dym (4.6) 

(27r)'^Vdet J_oo J-oo 

where Y = [yi, . . . ,ym)- We have the following extension of the formula (3.6). 
Proposition 4.1 

det - nr=i(i - tf)^ni<.<,<^(i - ut,r ^ 

The proof of Proposition 4.1 uses the following lemma. 

Lemma 4.2 Let fn{t) {n > 3) be any random polynomial and ti, . . . ,tm be any m 
real numbers. Let be the covariance matrix of the Gaussian random vector 

ifnih), /n(^l)) • • • ) fnitm), fnitm)) 

which is defined in (4-2) and (4-3). Then 

detA^) = P„(ti,...,t^) H {ti-tjf (4.8) 

l<i<j<m 

where Pn{ti, . . . ,tm) is a polynomial. 

Proof. To simplify notation we drop the indices m, n in the matrix A^ • We have 

A — {Aij)ij^i^___^rn 

where 



_ /E/,(t,)/n(t,) E^(t,)/;(t,)\ .4... 
~\^m)fn{tj) E/;(t,)/;(t,); ^^--^^^ 

In the following discussion we consider linear transformations of the matrix A which do 
not change its determinant. By substracting the first and second column of Aji from the 
first and second column of A^j, respectively, we get the matrix A^-*^) with the 2x2 blocks 

(1) _ /E^(tO(/n(t,)-/n(tl)) E^(t,)(/;(t,)-/;(tl))\ 

~ \E mKfnitj) - fnih)) ^ f'niti)U'n{h) - f'nitl)) ) ^ ^ 

Since fn{t) is a polynomial, we can take the factor [tj — ti) out of the first two columns of 
the matrix A^^\ and this proves that detA is divisible by {tj — tiY . Repeating the same 
operation on rows we get the factor [tj — t\Y. How to get {tj — ti)®? To that end let us 
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subtract the second column of the matrix A\1^ = An multiphed by (tj — ti) from the first 
column of the matrix A^^\ This produces the matrix 

(2)_ /E^(tO[/n(t,)-/n(ti)-(t,-ti)/;(ti)] E^(to[/;(t,)-/;(ti)]\ 
~\^m)[fn{t,)-fn{ti)-{t,-h)ah)] Efi{u)[fi{t,)-f^{h)]) ^^-''^ 

Now we can take {tj — ti)^ out of the first column and {tj — ti) out of the second column 
of the matrix A*^^-*. Repeating the same operations over the rows we get that det A is 
divisible by {tj — ti)^ . Finally, let us observe that by the Taylor formula 

fn{t,) - fn{tl) - {tj - = + 0{\t, - tlf), 

and 

fnitj) - fLitl) = {tj - tl)f:{t,) + 0{\tj - 

hence if we subtract the second column of the matrix A^^-* multiplied by {tj —ti)/2 from 
its first column, the difference is of the order of \tj — ti\^, and we can take the factor 
{tj — ti)^ out of the first column and {tj — ti) out of the second column. This gives the 
factor {tj — ti)^. The same factor is taken out of the raws, hence det A is divisible by 
{tj — ti)^. Similarly, it is divisible by {tj — ti)^ for all i ^ j, and hence it is divisible by 
their product. Lemma 4.2 is proved. 

Proof of Proposition ^.1. By (4.1), we have 

det A^ = lim detA^;^^ (4.13) 

n— >oo 

for all ti,...,tm in the interval (—1,1). In fact, the limit (4.13) holds for all complex 
ti, . . . , in the unit disk {\t\ < 1}, and the convergence is uniform on every disk {\t\ < r} 
where r < 1. Hence by Lemma 4.2, 

det Am = H{tu...,tm) n (^^-^j)^ (4-14) 

l<i<j<m 

where H{ti, . . . ,tm) is holomorphic in the unit disk. 

Now, let us consider the expression of det A^ in terms of the matrix elements 5ij of 
A^, that is 

det Am = ^e^Si^(^l) .. .Sna{2m) (4-15) 

where a is a permutation of {1, . . . , 2m} and = ±1 depending on whether a is even or 
odd. The common denominator of the sum in (4.15) is 

m 

11(1-^')' n (1-^^^.)' (4-16) 
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Therefore by (4.14), 



det Arr,, = 



nZii^-tirni<i<j<mi^-ut,y 



(4.17) 



where C(ti, . . . , tm) is a polynomial of ti, . . . , t^. Observe that (4.17) holds for all points 
ti, ■ ■ ■ ,tm in the unit disk, and so it can be extended to the whole complex plane. We are 
going to show that C{ti, . . . ,tm) is a constant, and moreover, that 



C{ti,. . . , tm) — 1 



(4.18) 



Let us look at the asymptotic behavior of det as ti ^ cxd while t2, ■ ■ ■ ,tm are fixed. To 
see it more clearly, let us change A^, to A^"* by substracting the {2i — l)th column and 
row from (2z)th column and row, respectively, for i = 1, . . . , m. Then 

\ 



A^i) = 




^(1) 

^12 
^22 


. AS 
. A« 






^m2 • 


'-^mra 



(4.19) 



where 



for A; = 1, . . . , m, and 



^kk 



i-ti 





1 



2^ 



(4.20) 



aW = 



2titj+tjt2_2t|-2f2 + l 



1 t-itj 

ti-ti 



(4.21) 



for z 7^ j. The leading powers of the elements of A^'', as ti — > oo, are 



A(2) = 



/l/tj 1/ti l/h 

1/tf 1/t? 1/t? 

1/tl l/tf * * 

1/tl 1/tl * * 



1/tl 1/tl * 



* 
* 



1/tl l/tl\ 

1/tl yti 
* * 



\i/ti 1/tl * 

where *'s stand for the terms of the order of 0(1). Therefore 



* 
* 



(4.22) 



det A, 



O 



ti 



oo. 



(4.23) 
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By (4.17), 



det A, 



const. C(ti, . . . ,tm) 



ti — > oo, 



hence C(ti, . . . ,tm) is constant in ti. The same argument on t2, ■ ■ ■ ,tm shows that it is 
independent of t2, ■ ■ ■ ,tm, so it is indeed a constant, say Cm, i-e.. 



det 



^rnY[l<i<cj<mi^'>' 



(4.24) 



To prove that Cm = 1, let us consider the asymptotic behavior of det A^ as ti — 1 with 
t2, ■ ■ ■ ,tm fixed and close to zero. Then on the one hand, we have from (4.28) that 



lim(l-t2)4detA^ = -^ 



det A 



On the other hand. 



A(i) 







m—li 



* 



(4.25) 



v * 



Am-1 



(4.26) 



/ 



where the terms * are regular at ti = 1. Hence the leading term of the Laurent series of 
det A^ at = 1 is (1 — ti)~^ det A^_i, which shows that 



Cn 



Cm-1 



det A^_i = det A^_i 



(4.27) 



Thus Cm = Cm-1- Repeating this argument we get that 



C'm — Cm—1 



C2 = Ci = 1 



Therefore Cm = 1- Proposition 4.1 is proved. 

Similarly we prove the following proposition. 
Proposion 4.3. 



( ^11 ^12 
0^21 "^22 



Vc 



with 



COri. — 
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(4.28) 



(4.29) 



and 



Put now 



Then by (4.29) and (4.30), 



ti = tanh(7rsi) i = 1, . . . ,m. 



iO. 



• = (1 - ^2)3/2(1 _ ^2^3/2 V^^^i^i ^"^^^ ' ^fc) Hfc^, COth^ 7r(g,- - gfe) 

' ^ cosh7r(si - Sj) 



In addition, by (4.7), 



det = 



ni<i<j<„^tanh^7r(si-Sj) 



a=i(i-^f)^ 



By (2.14), 



A- Csi s ^ - -^"^(^1' ■ ■ ■ ,^m) / V 



7r(l-t2) 



Now let us substitute y^'s in (4.6) by 



(1 - t2)3/2 



i = 1, . . . , m; 



then 



Km(tl: • • • ; ^m) — 



X 



2-ni<i<i<n,tanh"7r(si-s,) 

/oo /"CX) 
-oo J —oo 



where 



f C^ll ('■12 
C21 C22 



C2m 



with 



_ rifc^i coth^ 7r(gi - gfc) Hfe^j coth^ 7r(gj - gfc) 
cosh7r(si — Sj) 

Now substitute j/^'s in the formula (4.36) by 

Vi 



rij^i coth^7r(si-Sj) 
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i = 1, . . . , m; 



then 



J —OO J —OO 



where 



Thus 



F — 
^ m. — 



/OO fOO 
■ ■ ■ I \yi . . .ym\e~^^^^"''^^ dyi . . .dyra- 
-OO J — OO 



(4.40) 



(4.41) 



(4.42) 



Theorem 1.2 is proved. 

Proof of Theorem 1.3. Let 

/n(i) = Co+Cit+--- + C„_ir-^+Cnr, gn{t) = Cof + C^f~^ + ■ ■ ■ + Cn-lt + Cn. (4.43) 

Then if ^ is a zero of /n(t) then t^^ is a zero of gn{t). Hence if 1 < a < 6 then 

e/(a,6) = e,(a-\6-i), (4.44) 

where $,f{a,b) is the number of zeros of fn{t) in the interval (a, 6). Observe that the 
distribution of the vector of coefficients (cq, . . . ,c„) coincides with the one of the vector 
(cn, . . . , Co). Hence 



E [0-(ai, 6i) . . . ^f{am, 6m)] = E [^,(ar\ &r') • • • ^sfc , O] 

= E[eHar\6r^)...0(am\O]- 
Take = tj and 6j = tj + Atj, j = 1, . . . ,m, and get Atj — > 0. Since 

\a-^ -b-^\= a-> - 6| + 0(|a - bf), a ^ b, 



(4.45) 



we deduce that 



Hence 



Knmjti ■>■■■■) tm ) 
C]_ . . . 



Knmifl: • • • ; ^m)- 



because 



Knmi^X 7 • • • 5 ) -^nm(^l7 • • • 5 ^m) 

p(tr^)---K^m ) ~ P(^l)---P(^m) ' 

Pit-') = ,i_\-2| = 



(4.46) 



(4.47) 
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This proves that 



"'nm^^l) • • • ) 



^m) — "'nml.Sl) • • • ) 




Theorem 1.3 is proved. 

Proof of Theorem 1.4- Let . . . , < 1 and . . . , > 1. Denote by 

Dnmixi,yi, ■ ■ ■ , Xm, Urri'i ti, ■ ■ ■ tm) the joint distribution density of the vector 

ifnitl),fnitl), . . .,fnitl),fnitl),9nitl+l),g'nitl+l), ■ ■ ■ Qnitm) , Qnitm)) ■ 

Then 




(4.48) 



X -Dnm(0, yi, . . . , 0, ^1, • • • , • • • , tj')dyi ■ ■ ■ dy^. 



The covariance 



^ fn{ti)gn{t^^) = + '^7' + • • • + Ut]-'' + tj", 



1 < i < i < J < m, 



goes to as n — > oo, together with the partial derivatives in tj, tj, while 



lim Efn{ti)fn{tj) 



1 



l<i,j<l, 



1 titj 



(4.49) 



n 



lim Egr,{t-^)gr,{t-^) 




1 ' 



I + 1 < i,j < m. 



This proves that the limiting Gaussian kernel Dm = 



lim^^oo Dnm is factored, 



Dm{xi,yi, . . . Xm, ym;ti,---, ti,ti+i, ■ • • , ) = A (/) (a^i , 2/1, ■ ■ ■ , , |/z ; ti , . . . , ti) 

X Dm-l{g){xi+i,yi+i, . ..,Xm,ym\ti+i, ■ ■ -^tm)- 



Hence the correlation function Km is also factored. Theorem 1.4 is proved. 
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5. Correlations between Zeros of the SO (2) Random Polynomial 



In the following two sections we discuss the correlation functions of real zeros and 
the variance of the number of real zeros of the SO (2) random polynomial. Let fn{t) be a 
SO (2) random polynomial, that is 

n 

Ut) = Y.Ckt\ (5.1) 

fc=o 

where {c^} are real independent Gaussian random variables with 

Ecfc = 0, Ecl^al^ll). (5.2) 



In this case (1.4) reduces to 



An{t) = {i+er, 

B^{t) = nt{l + fT-\ (5.3) 

Cn{t) = n{l + ne){l + eY~\ 

which gives 

(see [BS], [EK], [BBLl], and others). The average value of real zeros is 

/oo 
Pn{t)dt = ^/^, (5.5) 
-oo 

and the normalized density, 



Pn{t) 



(5.6) 



^^ 7r(l + t2)' 

is the Cauchy distribution density. Observe that both (5.4) and (5.5) are exact relations 
for all n. Let us compute the two-point correlation function i^n2(^i,^2)- 

Define 

Then the real zeros of (fin{t) coincide with those of fn{t)- Hence, similar to (2.6), we can 
write that 

f'CO f'OO 

Kn2{ti,t2)=l I \yiy2\Dn2{0,yi,0,y2;ti,t2)dyidy2, (5.8) 
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— oo •/ — oo 



where Dn2ixi,yi, X2, j/2; ^15^2) is the distribution density function of the vector 

{(Pnih), ^Pnitl), ^n(^2), ^n(^2))- 



Observe that 



and 



^ 'fin{tl)^n{t2) ^ Pn{tl,t2) 



{i + ht2Y 



(l + t2)n/2(i+^2)n/2' 
{t2-tl) 



(5.9) 



E^;(t,Vn(t2) = ^(tl,t2) =np. (l+M2)(l+tD ' 
d^Pr 



-ri^Pn 



dtidt2 

{t2-h)' 



(5.10) 



(l + M2)2(l + t?)(l + ti) 

Define the random variable ipn{t) as 



+ npr 



{i + ht2)^ 



(5.11) 



Let Dn2{xi,yi,X2,y2',ti,i2) be the distribution density function of the random vector 
i^niti): ^niti), Vn{t2) i ^nih)) ■ Then after a change of variables, (5.8) is rewritten as 



Kn2{tl, ^2) _ ^2 
Pn{tl)Pn{t2) 



00 POO 



\yiy2\Dn2{o,yi,o,y2;ti,t2)dyidy2. 



(5.12) 



00 •/ —00 



From (5.10), 



EllJnitl)(finit2) = Pn 



EV'n(*l)V'n(^2) = Pr 



Vn{t2 - ti) _ 

(l+tlt2) ' 

n(t2-ti)' ^ (l+tDil + tl) 



(5.13) 



(l+tlt2)2 (l+tlt2)2 

hence the covariance matrix of the vector (<^n(^i)) V'n(^i)) ¥'n(^2)5 V'n(^2)) is 
/ 



1 




Pn 

\Pn (i+titj) ^" 



v^(t2-ti) 
(l+tit2) ^" 

^'^ (l+tlt2) ^ ^ 





1 



(5.14) 



where 



, ,^ n(t2-ti)^ (l+tf)(l+ti) 

" = "^'^''^^ = "(l + tlt2)^+ (l+tlt2)^ 
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(5.15) 



Suppose that ti and t2 are two distinct fixed points. Then as n — > oo, the quantity 

(5.16) 



\Pn\ = 



goes to exponentially fast, and hence A„ approaches the unit matrix exponentially fast. 
By (5.12) this implies that 

iim — - — - — - — - = 1, 

n^oo Pn{tl)Pn{t2) 

and the rate of convergence is exponential. In the same way we obtain that if ti, . . . ,tm 
are m distinct fixed points then 

Knm{ti, ■ ■ ■ ,tm) _ f _f , _ y^\Dn^{0, y^, . . . 0, y^;ti, . . . ,tm) dyi . . . dym, 



— oo •/ — oo 



Pnih) ■ ■ -Pnitm) 

(5.17) 

where Dnm{xi, yi, ■ ■ ■ Xm, ym'i ti, ■ ■ ■ , tm) is a (2m) x (2m) Gaussian density with the co- 
variance matrix 

= {^ij )i,j = l,...,mj 

where 

/i n\ / 1 ^Mti-tJ)\ 

aL"'= J J), '^'=^-'('"*MvEi^ ("8) 

where Pn(^i) ^2) and a{ti,t2) are defined in (5.9) and (5.15), respectively. For fixed different 
ti, . . . ,tm the matrix A„ approaches the unit matrix exponentially fast, and this implies 
that 

lim — - — — - = 1, 

n^oo pn{ti) . . .Pn{tm.) 

and the rate of convergence is exponential. This means the independence of the distribution 
of real zeros at distinct fixed points. 

The formula (5.17) is simplified if we make the change of variables t = tan 9 (stereo- 
graphic projection). Consider therefore the random function 

n 

9n{0) = J2ckSm\e)cos^-\e), -| < ^ < I . (5.18) 

i=o 

Then 

^n(^) = cos"^/„(tan^), (5.19) 
hence if {tj} are zeros of fn{t) then 

{r]j = arctauTj} (5.20) 
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are zeros of gn{0). When = 0, gn{0) has an extra zero r] = 7r/2. Since the probabihty 
that c„ = is equal to zero, the joint probabihty distributions of zeros of the functions 
gni^) and /n(tan6') coincide. By (5.6) the zeros {rjj} are uniformly distributed on the 
interval [— 7r/2, 7r/2]. The function ^n(^) is a Gaussian random function with zero mean 
and the covariance function 

Egn{ei)gn{e2) sin'^ Oi cos"-'^ Oi sin^ 02 cos"-'^ O2 = cos"(^i - 82). (5.21) 

The function gn{d) is periodic of period tt if n is even, and it is periodic of period 2tt if n is 
odd. It is convenient to consider gn{0) on the circle = R/(27r)Z of the length 27r. This 
circle is the covering space of the original circle R/ttZ. If gn{0) = then g^(0 + n) — as 
well. On the distribution of gn{0) is invariant with respect to the shift 

9 ^ a + 9 mod 2it, 

hence it is S0(2)-invariant. 

Let Knm{9i, . . . ,9m) be the correlation function of the zeros {rjj} of gn{9). Assume 

that 

9j -9k^0 mod TT, l<j,k< m. (5.22) 

Then (5.17) gives that 

Knm{9i, . . . , 9m) 



••■ / \yi . . .ym\Dnm{^,yi, ■ ■ ■^,ym;9i, . . . ,9m) dyi . . .dym, 

-00 J — 00 

(5.23) 

where Dnmixi, yi, ■ ■ ■ Xm, y-m] 9i, . . . , 9m) is a (2m) x (2m) Gaussian density with the co- 
variance matrix 

— {^ij^)i,j = l,...,m, 

where 

AH = cos"(^,-^,)f I , 2,f'^T,^'~H^a a^) (5-24) 

y ^ 3J \ y/Eta.n{9j - 9i) -niaja^{9i - 9j) + cos~^{9i - 9j) J ^ ' 

Observe that Dnmixi,yi, ■ ■ ■ Xm, ym',^i, • ■ ■ ,^m) is the probability distribution density of 
the vector 

n (9 ^ ^"^^^^ n (9 ^ ^^^^ 
gny^l), , • • ■■,gn\^m)i 



'n vn 
and it is nondegenerate provided that n > 2m — 1 and (5.22) holds. 

Consider now the scaling limit of the correlation functions. The straightened zeros 

are 

TT 
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They are uniformly distributed on the circle of the length 2^Jn. The limit m-point corre- 
lation function of {Cj} is 



kjn^sx,. . .,Sm) = lim 



-^nm(^l7 • • • 7 ^m) 



Let us find k^isi, ■ ■ ■ , Sm)- We have that 



(5.25) 



lim cos ' 



(5.26) 



and by (5.24), 



with 



lim = A = (Ay)jj=i,...,, 



1 7r(sj - Sj) 

7r{sj - Si) 1 - 7r^(sj - Sj)^ 



(5.27) 



By (5.23) this gives that 
k. 



/OO nOO 
■■■ \yi ■ ■ ■ym\dm{0,yi, ■ ■ ■0,ym; Si, . . . ,Sm) dyi . . .dym, 
-OO J — OO 

(5.28) 

and dm{xi,yi, . . . Xm, ym', si, ■ ■ • , 'SJJ^) is a Gaussian density with the covariance matrix A. 
Let Qhe mxm matrix obtained by deleting all odd rows and odd columns from the matrix 
A~-^. Then we can write (5.28) as 

km{si, ...,Sm)^ 2"^(det A)i/2 J "J \yi---ym\e~^^^^'^^dyi...dym. (5.29) 



In Appendix C below we prove that 

A > if Siy^Sj,l<i<j<m, 

so the formula (5.29) is well-defined when the points Si are distinct. 
For m = 2, (5.29) reduces to 



(5.30) 



2 pOO pOO 

^tA)v^ y_oo j- 



k2{si,S2) = 



47r2(det 



\y1y2\e 2 



dyidy2 



(5.31) 



where 



A = 



/ 1 





2 





1 -se-^'/2 (1 - s2)e-«V2 

e -se-'^/\ 1 



\se-'/^ (I_s2)g-.V2 








1 



s = 7r(si-S2), (5.32) 
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and 

17=(A-1){2,4}, 

i.e., O is the 2 x 2-submatrix of at the second and the fourth rows and columns. 
Observe that 

det il = — - — = — - — - — . 
det A det A 

A direct computation gives 

det A = (1 - e-*' f - s^e-^' (5.33) 



and 



where 



Since 



Q = 



A B 
B A 



^" det A ' det A " ^^"^^^ 



/oo poo zL / A \ 

j |yiy2|e-^(^^'+^^^'+'^^^^^) dy^dy2 = + -^f^ arcsinjj , 

where 5 = B/A (see Appendix A), we obtain from (5.25) that 

= (detA)VMetO + """^ 

(detA)V2 f 5 \ , , 

^ ' ' 1 + ^=^= arcsm5 (5.35) 



[(1 _e--')2-s4e-^Y/2 / 5 

' ^ ' 1 + , = arcsm 5 



where s = 7r(si — S2) and 

e-^'/^(e-^'+g^-l) 

1 — s2 2 — s2 

1 — e * — s e * 



(5.36) 



It is worth to note that Hannay [Han] has calculated the limit two-point correlation function 
of zeros for the complex random SU(2) polynomial, and our calculation of (5.35) is very 
similar to the one of Hannay. 



As s ^ 0, 



8 2 

det A = ^ + O(.^O), = 1 - ^ + 0{s'), 
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which imphes that 



\Sl - S2\ , _ _ |2 



+ 0{\si-S2n, Si -52^0. 



(5.37) 



As s — > oo, 

det A = 1 - s^e-'" + 0{e-'"), 5 = -s^e-^'/^ + 0{e-'"/''), 
which implies that 



k2{si, S2) = 1 + 



+ 0{{si - S2fe-^ (^1-^^) ), |si - S2H 00. 

(5.38) 



Thus we have proved the following theorem. 

Theorem 5.1. Let {tj} be zeros of a random S0(2) polynomial fn{t) of degree n, 

^Jn arctanTj 



and 



TT 



he the straightened zeros. Then the limit m-point correlation function km{si, . . . , Sm) of 
{Cj} given by the formula (5.29) where /S. is a (2m) x (2m) symmetric matrix which 
consists 0/ 2 X 2 blocks defined in (5.26), and O = (A~^)even is the my. m matrix of 
the elements of with even indices. The 2-point correlation function is given by the 
formula (5.35), and its asymptotics as si — S2 — >■ and |si — S2I — >■ 00 are given in (5.37) 
and (5.38), respectively. 

The graph of A;2(0, s) is shown in Fig. 2. 




0.2 0.4 0.6 0.8 1 



Fig 2: The two-point correlation function of straightened zeros of the SO(2) random polynomial. 
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6. Variance of the number of real zeros of the SO (2) random polynomial 



Here we calculate the variance of the random variable ^n(a, b) in the case when fnit) 
is the SO (2) random polynomial. By definition, 

Var Cn{a,b) = EC{a,b) - iEU(^,b)f. 



By (2.17), 



Since 



we obtain that 



EC{a,b)= / Pn{t)dt+ / / Kn2{ti,t2)dhdt2. 
J a J a J a 

ph ph 

(Een(a,6))2= / Pn{tl)dti / Pn{t2)dt2, 
J a J a 



(6.1) 
(6.2) 
(6.3) 



pb pb pb 

Vai^n{a,b) = / Pn{t)dt+ / / [Kn2{tl,t2) - Pn{tl)Pn{t2)]dtidt2 
J a J a J a 



pb pb pb 

= / Pn{t)dt+ / / 
J a J a J a 



Kn2{tl, t2] 



(6.4) 



Pn{tl)Pn{t2)dtidt2 



.Pn{tl)pnit2) 

When ti and t2 are separated, the difference 

Kn2 {ti,t2) _ ^ 

Pn{ti)pn{t2) 

is exponentially small, hence the main contribution to the last integral comes from close 
ti,t2- Let us put 



h = t, t2 = t + 



Pn{t)' 



Then 



r Kn2{ti,t2) _ ^ 

J a J a [Pn{tl)pn{t2) 



hence 



Thus 



where 



Pn{tl)Pn{t2)dtidt2 ^ / Pn{t)dtl (A;2 (0, s) - l)(is, 

J a J —oo 

pb r poo 

Varen(a,fe)~ / Pn{t)dt 1- / {I - k2{^,s))ds 

J a L J —oo 

Var ^„(a, b) - C^n / , 



/oo 
(A;2(0, s)-l)ds, 
-oo 

and /c2(-si, S2) is the two-point correlation fumction given in (5.35). In particular, 

Var ^„(— 00, 00) ~ C^/n. 

Numerical value of C is 



C = 0.5717310486902.... 
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Appendix A. Calculation of an Integral 

In this appendix we show that 



/OO /"OO 4 / A \ 

J_ \ym\e-'^^^y"+''y"-+^^y^y-Uyidy2 = (^1 + ;^^=^ arcsinjj , (A.l) 

where 6 — B/\/AC. By a change of variables we can reduce the integral (A.l) to 

/OO /-cx) 
-OO J —OO 



Changing then 



1 / 

y2 = ~^2), 



we obtain that 



-| poo /*oo 
J —OO J — OO 



Let 



a;i = , cos 9, 
^/l + 5 
r 

X2 = , sin^. 



Then 



2 roo 2 /•27r 



(1 - ^2)=^/^ Jo 
Evaluating the last integral we obtain that 

/= 1 

(1 - 52)3/2 

hence (A.l) is proved. 



cos^ 9 siv? 9 
2(1 - 52)1/2 • - - YTI ~ 

1 f^^ 

--^ / I (1 - 5) cos^ 9 -{1 + 5) sin^ 9\d9 

(1 - d^)'^/^ Jo 

1 f"^^ 

'5 -cos 29\d9 



d9 



(\/l — 5^ + 6 arcsin5) 
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Appendix B. Positivity of the Covariance Matrix 



Assume that fnit) = cq + cit + • • • + Cnf^ is a polynomial with independent random 
coefficients such that Ec/. = and < Ec^ < oo. We show in this appendix that the 
covariance matrix of the random vector 



Fn — {fn{tl)i /n(^l)) • • • ) fn{tm)^f'n{tm)) 



(B.l) 



is positive, if ti, . . . ,tm are pairwise different and n > 2m — 1. Consider the real-valued 
quadratic form associated with the covariance matrix, 

Q{a,P) = J2[»J<^kEfn{tj)fn{tk) + ajPkEfn{tj)f^{tk) 



= E 



k=0 



The generalized Vandermonde matrix 
is nondegenerate, hence 



,m; fe=0,...,2m— 1 



fc=0 



n 2 



3 = 1 



>0, 



provided that not all aj,Pj are zero. This proves that the quadratic form Q{a,P) is 
positive, hence the covariance matrix of the random vector (B.l) is positive as well. 

The proof remains valid if Ck are, in general, dependent random variables with positive 
covariance matrix {Vki — E CfcQ)fc,i=o,...,n- In this case 

n m 

Q{a, P)=J2 ^kidkdi > 0, dk = Y, + /^i^"') > 

k,l=0 3 = 1 



provided that n > 2m — 1 and not all ctj , (5j are zero. 
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Appendix C. Proof of the Inequality A > 



Let A = iAjk)j,k=i,...,m where 

^ \{Sk-Sj) l-{Sj-Sk) 

We prove in this appendix that A > 0. Consider the complex- valued quadratic form 
associated with A, 

Q{a,/3) = ^e-(*^-*'=)'/2[a,afc + (s, - Sfc)^,^^ + {sk - Sj)/3jak + [1 - {sj - Skf]/3jPk]- 
Using the formulae 

V2n J-oo 
1 

{sj - Sfe)e-(^^-^^^'/2 ^ / e'^'^-"'>{-ix)e-'''/^dx, 

V27r J-oo 
1 

V27r J-oo 

we can rewrite Q as 

Q{a, (3) = -j= / dxe-""'/^ V [ajake'^'^-'"^'' + ajP^e'^'^-^^^^'i-ix) 
V27r J-oo - j^ 

+ Pjake'^'^-"'^''{ix) + Pjpke'^'^-"^^''x'^]. 

Define 

m m 

j=i i=i 

Then 

Q{a,p) = -= / \f{x)+ixg{x)\^e-^"/^dx. 

The function f{x) +ixg{x) is not identically zero provided that not all ocj, Pj are 0. Indeed, 
for every test function (fi{x), 

/oo 
{f{x) + ixg{x))(fi{x)dx = ^{aj(fi{sj) + Pj(fi'{sj)), 

where 0{s) is the Fourier transform of (fi{x). We can localize the function 0{s) near Sj and 
make the last sum nonzero. This proves that f{x) + ixg{x) is not identically zero, and 
hence Q{a, P) > 0. Hence A > 0. 
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